FREE MONOTONE TRANSPORT 



A. GUIONNETt AND D. SHLYAKHTENKO* 

Abstract. By solving a free analog of the Monge-Ampere equation, we prove a non- 
^sj ' commutative analog of Brenier's monotone transport theorem: if an n-tuple of self-adjoint 

non-commutative random variables Zi, . . . , Z„ satisfies a regularity condition (its conjugate 
variables , . . . , ^„ should be analytic i\\ Zi, . . . , Z„ and should be close to Zj in a certain 
analytic norm), then there exist invertible non-commutative functions Fj of an n-tuple of 
^ ' semicircular variables ^i, . . . , Sn-, so that Zj = Fj{Si, . . . , Sn)- Moreover, Fj can be chosen 

to be monotone, in the sense that Fj — !^jg and g is a non-commutative function with a pos- 
itive definite Hessian. In particular, we can deduce that C*{Zi, . . . , Z„) = C*{Si, . . . , Sn) 
I and W*{Zi, . . . , Z„) = L{¥{n)). Thus our condition is a useful way to recognize when an 

' n-tuple of operators generate a free group factor. We obtain as a consequence that the 

g-deformed free group factors rq(R") are isomorphic (for sufficiently small q, with bound 
depending on n) to free group factors. We also partially prove a conjecture of Voiculescu by 
showing that free Gibbs states which are small perturbations of a semicircle law generate 
free group factors. Lastly, we show that entrywise monotone transport maps for certain 
<^ ' Gibbs measure on matrices are well-approximated by the matricial transport maps given by 

■ free monotone transport. 
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1.1. On a notion of density in free probability. Let Xi,...,X„ be classical random 
CN ■ variables. Thus X = {Xi, . . . , Xn) can be viewed as a function defined on a measure space 
tJ- ! {Q,uj) with values in (say) M". Of special interest is the law of X, which is the measure on 

^ I obtained as the push-forward /i = X^cj. Very often, it is assumed that /i is Lebesgue 

absolutely continuous and the density p = dfi/{Y[ dxj) then plays a key role. For example, the 
density is involved in the classical definition of entropy (J plogp ]^(ixj), Fisher information 

X : (/ ^ n dxj) and so on. 

I In passing to the non- commutative case, one assumes that Xi,...,X„ are self-adjoint 

elements of some finite von Neumann algebra (M, r), where r : M — )• C is a normal faithful 
trace. In other words, Xj are self-adjoint operators on some Hilbert space H containing 
a vector ^ so that r(T) = (T^,^ satisfies r(P(X)Q(X)) = t{Q{X)P{X)) for any non- 
commutative polynomials P and Q evaluated at X = (Xi, . . . , X„). 

The only (partially) satisfactory extension of the notion of joint law of Xi,...,X„ to 
the non- commutative case uses the moment method. More precisely, one defines the non- 
commutative law of Xi,...,X„ as the functional fi which assigns to a non-commutative 
monomial P the value r(P(X)). 
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Unfortunately, no satisfactory replacement notion of density has been obtained so far with 
the exception of the case n = 1. In this case, spectral theory gives a suitable replacement; 
however, in absence of any commutation between Xi, . . . , Xn there is no satisfactory theory 
of "joint spectrum" even if X^'s all act on a finite-dimensional vector space. This causes 
a number of problems in free probability theory. For example, Voiculescu introduced two 
definitions of free entropy (x and x*)^ but for neither of them is there a simple formula 
analogous to the classical case as soon as n 7^ 1. 

1.2. Free Gibbs states and log-concave measures. One hint that gives hope that a 
satisfactory replacement for the notion of density can be found lies in the existence of free 
analogs of strictly log-concave measures. Log-concave measures on R" are probability laws 
having particularly nice density: it has the form exp{—V{xi, . . . , Xn))dxi • ■ • dxn, where V is 
a strictly convex function. 

It turns out that these laws have free probability analogs: If K is a "locally strictly convex 
non-commutative function" (see [GS09] for a precise definition or section 2.9; one example 
of such a function is V{Xi, . . . , Xn) = | ^ ^| + f3W{Xi, . . . , Xn) for a fixed polynomial W 
and sufficiently small f3), there exists a unique non-commutative law Ty which satisfies the 
"Schwinger-Dyson" equation [Gui06, GMS06] 

TviP^V) = Ty (g) Tv{Tr{j^P)). 

This law Ty is called the free Gibbs law with potential V. Here ^ and ^ are suitable 
non-commutative replacements for the gradient and Jacobian, respectively. 

Free Gibbs laws associated to such convex non-commutative functions are very well- 
behaved. It is worth mentioning that they arise as limits of random matrix models as- 
sociated to probability measures with log-concave density exp{—NTr{V{Ai,...,An))) on 
the spaces of N x N self-adjoint matrices. Alternatively, the law Ty can be characterized as 
the minimizer of relative free entropy 

Xv{rv) = Xi^v) -rv{V), 

where x{'''v) is Voiculcscu's (microstates) free entropy [Voi94, Voi02b]. 

Free Gibbs states have many nice properties with respect to free stochastic calculus, various 
free differential operators and so on (see e.g. [GMS06, GMS07, GS09]). Thus the fact that, 
without knowing what a non-commutative density is, we are able to single out a class of 
non-commutative laws which are similar in property to classical log-concave measures (which 
by definition are measures with a nice density!) strongly suggests that a non-commutative 
version of density could be found. Moreover, if it were to be found, such log-concave measures 
are natural first candidates for a detailed study. 

1.3. Monotone transportation maps as replacements for densities. There is an al- 
ternative way of talking about densities in the classical case, which is a consequence of the 

following theorem of Brenier: assume that fi is a measure on satisfying some technical 
conditions (Lebesgue absolutely continuous, finite second moment, etc.) Let u denote the 
standard Gaussian measure on M" . Then there exists a canonical monotone transport map 
0^ : IR" — )■ from v to fi. Here by a transport map from u to /iwe mean a map (f) satisfying 
(f)*!^ — //; and is called monotone if (f) — Vip for some convex function ip (roughly speaking, 
this means that the Jacobian of is positive-definite almost everywhere). Our point is that 
such a transport map carries all the information contained in the knowledge of the density 
d/i/Yldxj (indeed, this density is essentially the Jacobian of the map). A sign of the fact 
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that (f)^ is a good analytic object is that it continues to exist even when the density does 
not (for example, when fi is not Lebesgue absolutely continuous, but does not give mass to 
"small sets", cf. [Bre91]). 

It is thus very tempting to ask if there is any analog of Brenier's results in the case of non- 
commutative probability. Here one immediately faces a major obstacle: unlike in the classical 
case, there are many more non-commutative probability spaces than classical probability 
spaces. Indeed, up to isomorphism, the unit interval [0, 1] with Lebesgue measure is the 
unique non-atomic probability space. At the same time, there are many non-isomorphic non- 
commutative von Neumann algebras. In fact, there is not even a (separable) von Neumann 
algebra that contains all others [Oza04], thus there cannot exist a "master law" so that all 
laws can be obtained as "push-forwards" of that law (as is the case classically). However, as 
we show in the present paper, one may hope to obtain such a result in certain cases (and, 
indeed, the Brenier map only exists under certain analytic assumptions about the target 
measure, even in the classical case). 

1.4. The main results. In this paper, we give the first examples of existence of non- 
commutative monotone transport. We thus take the first steps in the study of this subject. 

To state our results, recall that the Gaussian law u can be characterized by the following 
integration by parts formula: for any / : M" ^ M", denoting by x the vector {xi, . . . , x„). 



where Jf denotes the Jacobian of /. In other words, x = J* I, where J* denotes the adjoint 
of J viewed as an unbounded operator L^(M", v) — )■ L^(M„xn, and / denotes the nx n 
identity matrix. 

Recall that a free semicircular family S = {Si, . . . , Sn) is a free-probability analog of an 
iid n-tuple of Gaussian random variables. The semicircle law can also be characterized by 
the formula of the form 



where ^* is the adjoint of the map ^ , which is a suitable non-commutative replacement 
for the Jacobian. 

Our main theorem shows that monotone transport exists in the non-commutative case, 
provided that the law of an n-tuple Z = {Zi, . . . , Z„) is "close" to the semicircle law, in the 
sense that it "almost satisfies" the equation Z = ^*(/), i.e., Z is close to in a certain 

analytic function norm || • ||a (see §2.1). The conditions on Z automatically imply that the 
law of Z is a free "log-concave" Gibbs law. 

Theorem 1.1 (Existence of monotone transport). Let S — {Si, . . . , Sn) be a free semicircular 
family. If Zi, . . . , Zn E {M, t) are non- commutative random variables such that there exists 
F = F* in the closure of A for the norm \\.\\a so that J*{1) = ^F is so that \\Z - J^*{I)\\a 
is sufficiently small, then there exist an Yi, . . . ,Yn E W*{Si, . . . , Sn) so that: 

• the law of Y — {Yi, . . . , F„) is the same as the law of Z = {Zi, . . . , Zn), so that 
W*{Z) ^ W*{Y); 

• W*{S) = W*{Y) {= W*{Z)), so that Yj can be viewed as "non- commutative measur- 
able functions" of Si, . . . , Sn (we in fact show that these functions can be taken to be 
non- commutative power series); 

• Yj = ^jG for some G e W*{Si, . . . , Sn) and ^Y is positive-definite. 
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Here 3l and ^ are suitable non-commutative replacements for the gradient and Jaco- 
bian of a function; these differential operators play key roles in free probability theory 
[Voi98, Voi02b, Voi02a, Gui09, Gui06, AGZIO]). In particular, the last condition is a kind 
of convexity requirement on G (since = ^ QiG is a kind of Hessian of G). 

By analogy with Brenier's theorem, we call F = S^G "free monotone transport" from the 
semicircle law to the law of Z. We can view the function F as playing the role of density; 
the fact that we are able to construct F in the space of certain analytic maps is a reflection 
of the fact that the law of Zi, . . . , Z„ has a "nice density" relative to the free semicircle law. 

By Voiculescu's change of variables formula, one can claim that for laws satisfying the 
assumptions of our theorem, 

= T ®T ® Tr(log ^ F^) + universal constant. 

We thus see that free entropy can be expressed by a simple formula in terms of the monotone 
transport map, and is concave in this map. 

Our theorem yields a number of isomorphism results, since it gives a rather general condi- 
tion of when an n-tuple of operators generates a free group factor. In particular, the theorem 
applies to free log-concave Gibbs laws with potential of the form ^ ^ X| -|- l3W{Xi, . . . , X„) 
if /3 is sufficiently small: 

Theorem 1.2. Let W he a fixed self-adjoint polynomial in n variables and set 

1/5 = 1 YTj=i -^j + Let Tj3 be the free Gibbs state with potential Vp. Then for sufficiently 

small 13, W*{Tp) = L(F„). 

This theorem partially answers in the affirmative a conjecture of Voiculescu [Voi06, p. 240] 
(the full conjecture involves arbitrary values of (5). 
We also prove: 

Theorem 1.3. Let rg{W') be the von Neumann algebra generated by n q-deformed semicir- 
cular elements of Bozejko and Speicher. Then there are numbers qo{n) > so that for all 
|g|<goW, r,(R")^L(Fj. 

This settles (for small values of q) the 20-year old question of the isomorphism class 
of such g-deformed von Neumann algebras. Furthermore, since I/(F„) is a factor, is non- 
hyperfinite, strongly sohd [OPIO] (thus is solid and has no Cartan subagebras), has the 
Haagerup property and the complete metric approximation property, (again, for the small 
values of q) our result can be viewed as the culmination of a number of results about these 
algebras, including [BS91, BKS97, Nou04, SniOl, Sni04, Shl04, Ric05, Shl09, KNll, DablO, 
Avsll]. 

It would be very interesting to extend our main theorem (which is limited by its hypothesis 
to laws "close" to the semicircle law), and to study the analytical properties of the map 
(these would be non-commutative analogs of properties of the density associated to a classical 
law /i). For example, could be used to at least formulate multi- variables analogs of various 
regularity results in single- variable free probability. Similarly, Ff^ could play a prominent role 
in the study of multi- variable random matrix ensembles. 

1.5. Free monotone transport as a hmit of classical monotone transport. As we 

mentioned before, the free Gibbs law Ty can be obtained as a limit of classical Gibbs measures 
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in the following sense: for any non-commutative polynomial P, 

1 



Ty(P) = lim E (N) 



^Tr(P(Ai,...,/l„)) 



where //y^^ is a measure on n-tuples of self-adjoint N x N matrices given by 
dny^^ = exp{—NTrV{Ai, . . . , An)) x Lebesgue measure. 

Let W he a. fixed self-adjoint non-commutative polynomial and set V = + W. If we 

denote by /i^^^ the Gaussian measure, there exists a unique monotone transport map F^^^ 
which pushes forward /j,^^^ into //(T^. This map operates on matrices entrywise. On the 
other hand, let F be the free monotone transport map taking the semicircle law to the free 
Gibbs state Ty- Using functional calculus, F induces for each a map of n-tuples oi N x N 
matrices (a kind of "matricial" transport). In Theorem 4.7, we show that F^^^ —F — )■ in L^. 
In other words, the "entrywise" transport is well- approximated by "matricial" transport given 
by functional calculus on matrices. This is somewhat parallel to Biane's result [Bia97] which 
shows that in the Gaussian case, the "entrywise" Segal-Bargmann transform is asymptotically 
"matricial" and yields the limiting transform in free probability theory. 

1.6. The one-dimensional case. In the case that n = 1, so that we are dealing with a 
single random variable, non-commutativity goes away. In this section, we explain the proof 
of our main result in the case n = 1; indeed, many of the formulae that hold in the general 
multi-variable case have simple ad hoc explanations in the one- variable case. The arguments 
in this section are subsumed by the arguments for a general n; however, we feel that their 
inclusion is justified as they serve to clarify our proof of the general case. 

We thus start with a semicircular variable X viewed as the operator of multiplication by 
X on L^([— 2, 2],r]), where rj is the semicircle law: dr]{x) = X[-2,2]{x) ■ ^\/4 — x^dx. We write 
t(/i) = J h{x)dr]{x). Thus r] is the unique maximizer of the functional 

(1.1) x(/^) = jj log \s - t\dii{s)dii{t) - j \edii{t). 

We fix a function W which is analytic on the disk of radius A and let V{x) — ^x^ + W{x). 
We now consider measure rjv which is the unique maximizer of 



Xvil-i) = JJ log \s - t\dfi(s)dfi(t) - j V{t)di^{t). 

If we write Ty for the linear functional Tv{h) — J h{x)dr]v{x), then Ty is the free Gibbs state 
associated to V. 

Our goal is to construct a function F which is analytic on a disk of radius A' for some 
A < A' < A, so that Y — F{X) has law rfy when X has law rj. In other words. 



;i.2) I h{F{x))drj{x) = I h{y)dMy): 



or, equivalently, F^r] = rjy. 

By definition Xv{{x + 5f)^rjy) < Xvijiv) for all real numbers 5 which implies that the 
unique minimizer r]y satisfies the Schwinger-Dyson equation 

j V'{t)f{t)dMt) = II ^^^{^drjy{s)dvv{t). 
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If we write J f{s,t) = ^^^4^' W(a;) = V'(x) and identify L'^([-2,2]'^,driv(s)dr]v(t)) 
with L°°([— 2, 2], driv)^L°°{[—2, 2], drjv), then the Schwinger-Dyson equation reads 

(1.3) rv{f^V)^Tv^TvUf). 

In the case that V is strictly convex on a sufficiently large interval (e.g., if W is sufficiently 
small), the Schwinger-Dyson equation determines Ty (or, equivalently, rjv) uniquely [GS09]. 

Note that (1.3) is equivalent to saying that ^*(1 ® 1) = !^V. It is not hard to see that 
^*(1® 1) = 2 J -^z^^driviy) (here and for the remainder of the section, the improper integrals 
are taken in the sense of principal value). Thus (1.2) is equivalent to 

2 / dr]v{y) — V'(x) Vx e supp(r^y) . 

J x-y 

Replacing x by F{x) and y by F{y) and remembering that F^rj = rjv gives us the following 
equation for F: 

^ I F(x)- F(y) '^'^^y^ = V'{F{x)) = F{x) + W'{F{x)), Wx e [-2, 2] . 

To solve this equation, we try to find a solution of the form F(x) — x + f{x) with / small if 
W is. Then the equation becomes 

(1-^) 2 / -i-^-^L_*,(,) = . + /(.) + wv + m), 

where, as before, we use the notation ^f{x,y) — ^^^^E^j^- 

Using the fact that 77 is the semicircle law, so that 2 J ^^dr]{x) — x, we find that (1.4) is 
equivalent to 

To deal with the difference quotient ^f{x,y) = ^^^^^^^^^ ^ we add and subtract the term 
fix): 



fix) 



-drj{x) 



+ f{x) + W\x + f{x)). 



x-yl + / f{x,y) 
We can substitute (1.4) into the right hand side, to get 

= f\x) [x + f\x) + W\x + /(x))] + f{x) + W\x + f{x)) 

^{xf)\x) + \{f{xf)' +[W{x + f{x))t 

The resulting equation is equivalent to the previous one under the assumption that f'{x) ^ 
—1 almost surely. 

Since the right-hand side is a full derivative, it is natural to try to express the left-hand 
side as a full derivative, too. Simple algebra shows that 
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/ ^(/(^) - //(^^ y)) 1 + Jf^^^ My) =d.j \og{l + /f\x,y))d^{y). 
Substituting this into (1.5) and removing derivatives from both sides finally gives us 

2 j \og{l + / f{x,y))dii{y) = {xf){x) + ]^{f{x)f + W{x + f{x)) + const. 

Finally, we seek g so that f — g' solves the previous equation. This gives us the following 
equation for g: 

(1.6) xg'{x) = -W{x + g\x)) - \[g\xyf + 2 ^ log (^1 + ^^^^^^^^ dr]{y) + const. 

The constant can be fixed by requiring that both sides of the equation vanish at a; = 0. Note 
that the operator ^ g = xg' multiplies monomials of degree n by n. Let S be the inverse 
of (defined on polynomials with zero constant term),E(^)(a;) = Jq y~^g{y)dy, and let 
g — J^g so that g — "Eg. With this notation, our equation becomes 

gix)^-W{x+{Egy{x))-^m){x)y + 2 J log (^i + iMM^^^^ dr}{y) + const. 

We can now rewrite this equation as a fixed point equation 

g{x) = ^,^{g{x)). 

Note that ^'a;(0) = W{x). For W small, this equation can be solved by iteration in the 
space of analytic functions converging on a large enough disk. The essential fact here is 
that g ^ "if xid) is Lipschitz for a certain analytic norm for any x. This gives us a function 
g, which is sufficiently small. Retracing our steps (and noting that for sufficiently small 
W, f is small, so that |/'(a;)| < 1 almost surely), we get back the desired map F given by 
F = x + f = G'ifG = lx^ + g{x). 

It is rather straightforward that if we consider Vp — + Wfs{x) for some family VFg 
analytically depending on (3 and sufficiently small, then our transport map F — Fy^ also 
depends on /3 analytically. 

It is worth noting that for W sufficiently small, the map F{x) = x + f{x) has a positive 
derivative and is thus monotone. It is therefore the unique monotone map satisfying F^rj = 
Tjv- In particular, the map F gives optimal transport between r) and rjv- 

It was pointed out to us by Y. Dabrowski that equation (1.6) can be regarded as a free 
analog of the classical Monge-Ampere equation. To see this, put V{x) = and V = V + W, 
and consider the classical Monge-Ampere equation for the monotone transport map H from 
the Gibbs measure d/j, — e:Kp{—V{x))dx to dp, — exp{—V{x))dx. Let JH stand for the 
Jacobian (derivative) of H. Then the Monge-Ampere equation reads: 

det{JH{x)) = """"^^'-^^^-"^{^ = ew{V{H{^)) ' V{x)). 
exp{—V{H{x)) 

(We write det JH even though JH is a 1 x 1 matrix to make the analogy clearer) . Taking 
logarithm of both sides and using that log det ^4 = Tr (log A) if ^4 > gives 

Tr{\ogJH) = V{H{x)) - V{x). 
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This is quite reminiscent of (1.6), if we rewrite it by setting / = g', F{x) — x + f{x) and 

2 1 log(/F)rfr;(|/) = |i(F(a;))2 + iy(F(a;))|-ia;2 = y(F(a;))-V^(a:). 
Later in the paper, we also consider the n-variable version of (1.6), which has the form 



(1 ® r + r ® l)Trlog = y 
where is a certain symmetrization operator. 

1.7. Organization of paper. To prove our theorem, we consider in section §2 a certain 

space of analytic functions in several non-commuting variables. This is the space in which 
our construction of the monotone transport map takes place. We also discuss the various 
differential operators which are the non-commutative replacements for gradients and Jaco- 
bians. We finish the section with a kind of implicit function theorem for non-commutative 
analytic functions. 

The next section is devoted to the construction of a non-commutative monotone transport 
map from the semicircle law r to the law Ty satisfying the Schwinger-Dyson equation with 
V — I ^ X| -I- 1^. We assume that 1^ is a given non-commutative analytic function (which 
we assume to be sufficiently small in a certain norm). 

The next section, §4, is devoted to applications. Finally, the last section collects some 
open questions and problems. 

Acknowledgments. The authors are grateful to Yoann Dabrowski for many useful com- 
ments and discussions. We also wish to mention that the idea of looking for a free analog 
of optimal transport (although via some duality arguments) was considered some 10 years 
ago by Cedric Villani and the authors; although the precise connection between optimal and 
monotone transport is still missing for the moment, such considerations have been an inspi- 
ration for the present work. We also thank an anonymous referee for his numerous comments 
who helped us to improve greatly our article. 

2. Notation and an Implicit Function Theorem. 

2.1. Non-commutative polynomials, power series and norms || • ||^. We will denote 
by ^ = C{Xi,...,Xn) the algebra of non- commutative polynomials in n variables. Let 
G he the linear span of polynomials with zero constant term. 
Following [GMS07], we consider on this space the family of norms || ■ ||^, A > 1, defined 
as follows. For a monomial q and arbitrary P e C{Xi, . . . , let Ag(P) be the coefficient 
of q in the decomposition of P as monomials; thus P — Xq{P)q. Then we set 

g:degg>0 

We'll denote by ^^^^ (resp., ^''^^) the completion of ^ (resp., ^) with respect to the 
norm || • \\a. This is a Banach algebra, and can be viewed as the algebra of absolutely 
convergent power series with radius of convergence at least A. 

Note that the norm || • ||^ has the following property: whenever Ti, . . . ,T„ are elements 
of some Banach algebra Q and ||Tj||Q < A, there exists a contractive map ^^"^^ — >■ Q 
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sending Xj to Tj. Hence, in the following we will see variables as elements of the set of 
polynomials However, will still also denote sets of maps, namely the completion 

of polynomials in the norm ||.||^. 

is the opposite algebra, with multiphcation defined by a°^.b°^ = {ha)"^. It is equipped 
with the same norm 

2.2. The operators ^ , T., U, ^ , d and ^. Let us denote by .yV the linear operator on 
S2f = C{Xi, . . . , Xn) that multiplies a degree k monomial by k. Let us also denote by 

the inverse of ^ pre-composed by the projection U onto given by 

H: P^P-P(0,0, ...,0). 
If 5f e we write ^jg for the j-th cyclic derivative of g [Voi02a]. For a monomial q, 

Note that S^jg = ^jUg. We'll denote by ^g the cyclic gradient (viewed as a vector): 

% = {^jg)]=,. 

If / = (/i, . . . , /n) with fj e then we'll write ^/ e M„xn(=2^ ®^'^) for the matrix given 
by 

where dj is the j-th free difference quotient [Voi98] defined as the derivation from ^ to 
£^ ® satisfying djXi — 5j=i, that is if g is a monomial djq — 'Ylfq=AXjB -A® B. 

2.3. Notation: Tr, ^* and We will often assume that r is a positive trace and also that 
1 ® 1 G dome?*, when dj : L^(^/, r) — )■ L^(^/, r)®L^(^/, r) is viewed as a densely defined 
operator. Under this assumption, £^ ® s^°p belongs to the domain of d* [Voi98]. We set for 
q^Yl,o.i®h ^ ^ ® and g e , 

i 

For q e M^xnl-s/ ® ^"p), and g = G we'll write 

Tr{q) = X %i e ^ (8) =2/°^ 




We will also denote by ^ multiplication in Mnxni-^ 8) ^2/°^) and write g^h = ^gihi if 
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2.4. Non-commutative notions of transport. Let (Af , (p) be a von Neumann algebra 
with trace 0, and let Xi, . . . , X„ G M be self-adjoint elements generating M; thus M can be 
regarded as a completion of the algebra and the trace induces on ^ a trace tx (called 
the non-commutative law of Xi, . . . , We write M = W*{tx). Let Yi, . . . , y„ e be 
self-adjoint elements generating another von Neumann algebra with trace ijj; we denote 
by ty the corresponding linear functional on so that N = W*{ty)- 

Definition 2.1. By transport from tx to Ty we mean an n-tuple of self-adjoint elements 
Yi, . . . ,Yn G M having the same law as Yi, . . . , Y^. We call such transport monotone if 
Y — {Yj)j belongs to the L^-closure of the set 

{^g:ge A\g(X^, . . . , X„) and ^% > 0}. 

(as part of the definition of monotone transport, we are making the assumption that the 
variables Xi, . . . , X„ are making the assumption that Xi, . . . , X^ are algebraically free, and 
so ^ Qig is well-defined for any g in the algebra these variables generate). 

When M is abelian, an element of M is an (essentially bounded) function of Xi, . . . , X„ 
(and n = L) So in the case that M is abelian, our definition reduces to the statement 
that the random variables Yi, . . . ,Yn are expressed as (bounded measurable) functions Yj = 
fj{Xi, . . . ,Xn)- Since the law of Y is the same as the law of Y, the push-forward of the 
measure tx via the function / = {fj)j is exactly the measure Ty- So our notion of transport 
coincides with the classical one in the commutative case. 

Our definition of monotone transport is analogous to the classical case where one requires 
that Y belongs to the closure of the space of gradients of convex functions. However, note 
that a different notion of derivative is being used. Nonetheless, when n — 1 our requirement 
reduces to asking that Y — f{X) for some / : R — > R so that f — g' for some g and 
f{x)-f{y) ~> g ■ g r -g monotone and q convex. 

x—y — 1 1 J o 

Note that if a transport from tx to Ty exists, then there exists a trace-preserving embed- 
ding N = Wiy) into M = W\X) given by Yj ^ Yj. 

2.5. The norm || ■ Let us denote by || ■ \a®^a the projective tensor product norm 
on (g) (^(^))°^'. In other words 

II y^fli ® hi\A®^A = sup \\Pl(S^ai O bi)\\ 

PI 

where the supremum is taken over all maps PI valued in a Banach algebra so that P/(l ® 6) 
and PI {a ® 1) commute and have norm bounded by the norm of b and a respectively. In 
particular, taking PI to be given by the right and left multiplication respectively, we see 
that for all g e ^ (8) J^/"^ and g E £/, we have 

ilgtt^lU < lklU®..i||^IU- 

We will also denote by the same symbol the norm on A'f„xn(-2/ ® (or its completion) 

given by identifying that space with the Banach space of (left) multiplication operators on 
® £/°P)". II ■ 11^^^^ also extends to (^ ® ^''p)" by putting for F = (Fi, . . . , F„) e 
^ (8) ^°P)" 

II-P'IU^^A = max ||Fj||^0^A. 

l<i<n 
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In consequence of its definition, we see tfiat the maps 

are contractive for the norm || • ||yi<g),rA on M„xn(=2/ (8) =2/*^) and the norm || • \\a on =2/". 

2.6. Cyclically symmetric functions. 

Definition 2.2. We say that a polynomial g in ^ is cyclically symmetric if for any ii, . . . , ik, 
the coefficient of the monomial Xi-^ ■ ■ ■ X^^ in q is the same as the coefficient of the monomial 

Let ^ denote the operator on defined on monomials by 

1 ^ 

^ Xi^ ■ ■ ■ Xi^ — - ^ Xi^^^ • ■ ■ Xi^Xi^ ■ ■ ■ Xi^. 

P r=l 

Clearly, sends ^ to the space of cyclically symmetric operators. Note also that S^g = 
S^gjj^X = '^■^igXi. Moreover, one clearly has that \\^g\\A < \\9\\a- Finally, note that 
9{yg) = 9g. 

2.7. Implicit function theorem. We next establish a kind of implicit function theorem 
for elements of the space [s^'^^^'Y. For later purpose note that if F e F is infinitely 
differentiable in s^^'^"> for any A' < A. These results are most likely folklore, but we were 
not able to find a precise reference. The proofs follow very closely the classical proof of 
the implicit function theorem using a contraction mapping principle. For F = {FjYj^-^ G 
(s/^^^)^, we will use the notation = sup^ We have denoted by ^ the algebra 
C(Xi, . . . , Xn) of non-commutative polynomials in n variables. Let us write =2/ * =2/ for the 
algebra of non-commutative polynomials in 2n variables, C{Xi, . . . , Xn, Yi, . . . , y„). Also, 
remember that variables and maps are all seen as converging power series. 

Theorem 2.3 (Implicit Function Theorem). Assume that F — {Fj)'j^-^ where 

Fj{X,,...,Xn,Y,,...,Yn) e 

is a non- commutative power series in 2n variables, and assume that 

F(0,...,0,0,...,0) = 0. 

Let Q = ^2^(0, . . . , 0) (here ^2 refers to differentiation in the variables Yi, . . . , Yn, with 
Xj 's held fixed), and assume that Q is invertible. Then there exists an G (0, A) and a 
finite constant C so that for all a e (0,ao), any U e ^^^^ satisfying \\U\\a < a, there exists 
a unique V G (^s/*^"^))" satisfying \\V\\a < Ca which solves the equation 

F(U, V) = 0. 

Moreover, there exists a function f G {s?/^'^^)"' so that V = f{U). 
Proof. Put 

G{U,V)^V-Q-'i^F{U,V). 

Then, since Q is invertible, 

G{U, V) = V ^ F{U, V) = 0. 
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Moreover, for all U, V since F(0, 0) = 0, 

G{U, V) = R{U) + H{U, V) 

with R{U) = —Q^^i^ ^iF{0, . . . , 0)i^U and H{U, V) is a non- commutative power series in 
[(^ * ^)*^"^'''] for any A' < A whose decomposition in monomials has only monomials of 
degree greater or equal to two. Therefore, for any U, U' so that ||?7'||a < a and any 

y, V so that \\V\\a < b, \\V'\\a < b, with a,b < A' < A, we have 

(2.1) \\H{U,V)\\a < CiA'Ka + bf 

(2.2) \\H{U,V)-HiU',V')\\A < CiA')ib\\U-U'\\A + a\\V-V'\\A) 

with a bounded constant C{A') which is an increasing function of A' < A. We next fix U 
so that \\U\\a < a and construct a solution by putting Vk+i = G{U,Vk) and Vq = 0. By 
construction, we have by (2.1) 

ll^fe+ilU < \\RiU)\\A + \\H{U, Vk)\\A < ||Q-'#/ii^(0, . . . , 0)\U^Aa + C{A'){a + WYkUf 

which shows the existence of a constant C so that ||Vfe||A < Ca for all /c > provided a < Qq 
for some oq > 0. Now by (2.2), we get for all A; > 1, 

11^4+1 - Vkh < C{A')a\\Vk - T4_i|U 

and hence Vk is a Cauchy sequence for C{A')a < 1. Therefore it converges when a < 
Qq a C{A')~^. The limit satisfies ||K|U ^ F{U,V^) = 0, which proves the 

stated existence. Moreover, by induction we see that for all k, T4 — fk{U) for some map 
fk £ and therefore the limit for \\.\\a satisfies the same property since .c/^'^^ is stable 

under composition. Finally, this solution is unique by (2.2) which guarantees that two 
solutions K and VI would satisfy 

||K - K'lU < C{A')a\\V. - K'lU 
implying that K = K when C{A')a < 1. □ 



2.8. Inverses to non-commutative power series. The following corollary shows that an 
"absolutely convergent non-commutative power series" has an "absolutely convergent inverse": 

Corollary 2.4. Let A > A' and consider the equation Y ^ X -\- ^fp{X) with e (^g/^^^)" 
and \\X\\a < A'. Then for sufficiently small /3, 3G e (=<^(^))" so that X = G'(Yg). 

Proof. We may without loss of generahty assume that fis{0) = by replacing Y with Y — 
^/^(O). U ^ ^Y,V ^ pX are solution of 

F{U,V) = U-V-/3'f^{(3-'V) = 0. 

Note that F(0,0) = and Q = ^2F{0,0) = -{I + /3 / f{Q)) is invertible for sufficiently 
small 13. Moreover, \\U\\a < (3A' + (3'^\\f\\A = (3 A" is small when 13 is and since fp G {s^^'^^Y, 
we can apply the implicit function Theorem 2.3 to conclude that the solution exists and is 
unique in \\V\\a < CpA". □ 
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2.9. (c, M)-convexity. Let V e =2/^ be a non-commutative power series. Following 
[GS09], we say that V is (c, M)-convex if for any n-tuples of self-adjoint operators X = 
{Xi, . . . , Xn) and Y = (Yi, . . . , F„) acting on a Hilbert space H and satisfying supj{||Xj||, \\Yj\\} 
M, V{X) is self-adjoint and the operator 

(W(X) - W(F)).(X -Y)- c{X - Y).{X - Y) 

is non-negative. Here we use the notation X.Y = XjYj + YjXj. 

3. Construction of the free monotone transport map. 

3.1. Outline of the proof of the main result. To prove existence of the transport 
map, we start with an n-tuple of free semicircular elements We then seek 

to find some elements Yi, . . . ,Yn in the von Neumann algebra M generated by Xi, . . . , X„ 
so that Yi, . . . ,Yn have law Ty, which is the unique log-concave free Gibbs law satisfy- 
ing the Schwinger- Dyson equation with potential V. We assume that V has the form 
V = + that 1^ is a "non-commutative analytic function" of Xi,...,Xn 

which is small in || ■ ||^. 

To construct Y = (Yi, . . . , F„) we seek / = (/i, . . . , /„) so that Y = X + f (i.e., Yj = 
Xj + fj). Assuming we could find such Y with fj analytic and of small norm, the Schwinger- 
Dyson equation for Y becomes an equation on /, of the form 

{rrjf) = ^ + ^ + + 

(Here ^* is the adjoint to Since for the semicircle law, ^*(1) = X, this equation 

further simplifies as 

^* (tTJi) + ^ + ^^^^^^ + /) = 0. 

One is tempted to solve this equation by using the contraction mapping principle; however, 
one faces the obstacle of bounding the map /i->- ^*(^//(l + ^ f))- Fortunately, under 
the assumption that / = Qig and further assumptions on invertibility of 1 + one is able 
to rewrite the left hand side of the equation as the cyclic gradient of a certain expression in 
g. This gives us an equation for g. It turns out that then one can use a contraction mapping 
argument to prove the existence of a solution of this equation for g\ this is the content of 
Proposition 3.13. We then retrace our steps, going from g to f — S>g and then toY — X-\-f, 
proving the main result of this section. Theorem 3.15. 

Fix C > 4 and let W = W* G ^'"^ be a cyclically symmetric function belonging to the 
closure of ^ in the norm 11 • 11^. Let / e and let Yj = Xj + fj, so that 

Y = X + f. 

The main result of this section is that under the assumptions that W is sufficiently small 
(i.e., and \\^W\\a(^„a are small), Y satisfies the equation 

(3.1) d*y.{l ® 1) = Y, + %.(W^(Fi, . . . , Fn)) 

or, in short 

J*{1) = Y + {^W){Y) 
where we denote by 1 the identity element of M^xniM ® M°^). This will imply that Y has 
law TiY^x^+w-i giving us the desired transport map (see Theorem 3.16). 
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To avoid confusion, we will use the following convention: all differential operators 
d) which either have no indices, or have a numeric index, refer to differentiation with 
respect to Xi, . . . , X„. Operators that involve differentiation with respect to Yi, . . . , F„ will 
be labeled as , dy^ , etc. 

3.2. A change of variables formula. 

Lemma 3.1. Put M = W*{£/,t). Assume that Y is such that = {dxyijij e 

Mn{M^M°P) is bounded and invertible. Assume further that 1 ® 1 belongs to the domain of 
d* for all j. 
(i) Define 

n 

1=1 

where # denotes the multiplication 

(a (8) b)^(A 0B) ^ Aa® bB. 
Then dj = dy, and dy.Xj = {{^Yy%. 

(n)d*y^{l^l) = Eid*x,iUY)-%. 

(Hi) Assume in addition that Yj = S^jG for some G in the completion of with respect to 
II • \\a- Assume that G = G*. Let (a ® b)^ = 6* ® a*. Then {JY)\- = {/Y)ij and {/Y)-^ 
is self-adjoint in Mn{M®M"P). 

Proof Let Q = JY . (i) We verify that 4^ = Y.iQxM{Q-^)io = EiQkMQ-% = 
lfc=jl 1 

To see (ii), we compute 

{d*y.{l^l),X,,...Xi^) = {l^l,dy.{X,,...Xi^)) 

= J](i®i,ax,(x,,---x,j#(g-i),,) 

e 

I 

To verify (iii), we first claim that if F = ^G for some G (not necessarily self-adjoint), 
then iy^Y)ij = cr[{^Y)ji], where a{a (8> 6) = h®a. Indeed, let G = Xj^ ■ ■ ■ Xj^.; in this case 

diYj= CA®B+ B® CA. 

G=AXiBXjC G=AXjBXiC 

Reversing the role of i and j amounts to switching the two sums, i.e., applying a. 
We next consider G of the form 

G = Xjj • • • Xj^ + Xjj, • • • Xj^ = P + P*, 

where P = Xj^ . . . Xj^,. In that case 

Yj = P^ + A*B* 

P=AXjB 
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and (3.4) 

P=AXjB BA=RXiS 

Thus, if we denote by * the involution (a (g) 6)* = a* 6*, 

(3.2) {JY% = a{{JY),,) = (^F),, 
proving that \s self- adjoint. Moreover we find that 

(3.3) {{/Y),,)^ = ^{{/Y)l,) = a{{/Y),,) = {/Y),,. 
where the last two equahties come from (3.2). 

□ 

A direct consequence of Lemma 3.1, (ii) and (iii), is the following. 

Corollary 3.2. Assume that g e 1 ® 1 G dom9* for all j, and put G = ^J^Xj + g. 

Let Yj — Xj + fj so that Y — S>G. Assume that 1 + ^ f is invertihle. Then Equation (3.1) 
is equivalent to the equation 



(3.4) [l + Jf J = ^ + ^ + (^^^^^ + 

3.3. An equivalent form of Equation (3.4). 

Lemma 3.3. Letr be the semicircle law. Assume that the map ^ i— )■ (^/)#C+C invertihle 
on (=e/^^))", and that f — S>g for some g — g*. Let 

K^-J*o ^ - id. 

Then equation (3.4) is equivalent to 



(3.5) K{f)^^{W{X-rf))-r 



Proof. Using the formula ^ = 1 — and the fact that ^*(1 1) = X, we see that (3.4) 
is equivalent to 

J! 

We will now apply (1 + ^/)# to both sides of the equation (this map is, by assumption, 
invertihle, so the resulting equation is equivalent to the one in the previous line): 



+ / + (W)(X + /) = 0. 



^ ' / + (W)(X + /) + (//)#^*^ 



+ ufw + ufm^w){x + /) = 0. 



,2 



Using that ^ x - we get: 
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Thus we have: 

K{f) = {^W){X + /) + + /) 

We now note that, because of cychc symmetry and the fact that / = with g — g*, we 
may use (3.3) to deduce: 

^(w(x + /)) = (^w)(x + f) + {jf)4{^w){x + /) 

Hence (3.4) is equivalent to 



K{f)^^{W{X + f)) + 
as claimed. 



□ 



3.4. Some identities involving ^ and It will be easier to give proofs of several 
identities involving ^ and ^ if we use the following graphical notation. We will denote the 
monomial Xi^ ■ ■ ■ Xi^ e by 



We will also write 



. . . . Zj, 



for the element of all of whose entries are zero, except that the j-th one, which is equal 
to X,, • • • X- . We will also write 



h ■ 






■Ji 



for the element Xi^-- - Xi (g) Xj^ ■ ■ ■ Xj e ^/ ® ^c/"^, and 



. 21 . . . tp ^ 
Jq---Jl 



for the matrix in Mnxni-^ ® ■^°) whose only nonzero entry is in the j'-th row and k-th 
column, and is given by Xi-^ ■ ■ ■ Xi^ ® Xj-^ ■ ■ ■ Xj^. 



j . . k 

Jq---Jl 



simply denotes the case where the entry [jk) is equal to 1 (E) Xj-^ ■ ■ ■ Xj^. 

With these notations, the operators ^ and ^ have nice graphical interpretations: 



= E 
= E 



• 1q-\-l ■ ■ ■ Ipll ■ ■ ■ 1q-l 





■h-l A 
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• ^<Z2 + 1 • 




*«2 




'92-1 ■ 





where the indices ii, . . . ,ip are cychcally arranged around the sides of the box. 
Let us use the convention 



3q---3i 



Jq' Jl 



6k- 



. ii . . . ipii ■ ■ - ip 



3q'---3l3q---3l 



Lemma 3.4. Assume that 1®1 belongs to the domain of d* for all j . Let g e s^'^^^ and let 
f — Sig. Then for any m > —1 we have: 



m + 2 



i=i 



Proof Let us say that g = Xlp ipO^hl-.tp^ii ■■■^ip and / = ^g. Let us write = 
9*(l(g)l). 

Recall [Voi98, Prop. 4.3] that ii B ^ Xi^ . . .Xi^® Xj^ ...Xj^es^® s^°^, then, if we 
denote by m the map a®h ® abc, we have: 

d*(A®B) = (A(g)S)#^^--mo(l(g)T(g)l)o(9j(g)l)(A(8)S) 
-m o (1 (8) T (2) 1) o (1 (g) dj){A (8) B) 

= A^jB - ^ dj=itXi^ . . . Xi,_jr(Xi,_^^ . . . Xi^)Xj^ . . . Xj^ 
t 

s 

Note that ^ ^ j)m+i^^._ It foHo^g ^^^^ 

J] {m o (1 ® r ® 1) o (9, ® mufr^%i) 



\m+l] 



} 



o (1 ® r ® 1) o (^z ® mUf) 

r 

5] {m o (1 ® r ® 1) o (1 ® dimj^fr+%i) 
I 

J2i^Mm o (1 ® r ® 1) o (1 ® dim^fr+%)] 
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We will concentrate on the first term, the other term being similar. It is not hard to see 
that, by the derivation property di and the fact that ^ f#{^ f)"^~^^ (^Z)™"^^, 



r 

= o (1 ® r ® 1) {[{di ® f)kr]M{J fr^%i) 

r 



where we use the notation a <S> b <S> c^2P 'S> q — a <S> bp <S> qc. Hence (using the notation 
a (g) — ab), we find 



Pi 



E E ■ 

(m+2) .(m+2) q^^ri=l 
h '•••''Pm+2 



Pm+2 

E 



m+2 

n 



a 



(Pu) 



9m+27^rm+2 = l "=1 



E 

qi<t<ri 



^gi+l • • • H-l -(l) 



,•(1) 



,•(1) 



■ ■ ■ ^n+i 




.(m+2) 




'*i2m+2 •(m+2) 


.(m+2) 


^9m+2-l ■ 


■ ''rm+2 + 1 



#1 



X5(m+2) .(I)T(X(I) ...X(1) • X(2) 



V2-1 



X 



(m+2) 
'''Jm+2 + 1 



...X 



(m + 2) 
'"m+2-1' 



.)■ 



Thus 



m+2 



w= E E ^fi.=.n 



q; 



(Pu) 



Pl...Pm+2 „. „. 



U=l 



E 

iji<t<ri 





•(1) 




(1) tW 




■ 'r-1+1 



•(2) 



''92 



^92-1 



^(2) 
^2 + 1 



A2) 

Or 2 



.(m+2) 

*'?m+2 •(m+2) 

9m+2-l 



,-(m+2) •(I) 
Vm+2 + l*t-l 



'?1+1 



#1 



X (5 



(m+2) 
'•m+2 " 



't+i 



'92 + 1 



(m+2) 
^+2+l 



(m+2) 
''■m+2-1' 



.)■ 
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It follows that 

w = EE ''.<:>^. n E 



m+2 



91 : '1 —"Pn 

Pl...Pm+2 „ u=l qi<t<ri 





•(1) 




(1) 


V-1 ■ 





.(2) 


.(2) 


.(2) 


.(2) lr2 


'92-1 • 


■ ■ V2+I 



• (m+2) 


• (m+2) 






^9m+2 •(m+2) 


• (m+2) 'l'rm+2 


4" 




?m+2-l ■ 









#1 



Xr(X.(l) . . . X.n) • • • X.(2) . . . X.{2) ■ • ■ X.(m+2) . . . X.(m+2) ) 
«; I 1 « 1 111 t ^ 111 It' 

"t+l Vj-l '92 + 1 V2-I 9m+2+l '•m+2-1 



PI 



m + 2 



E E E 

Pl---Pm+2 jU) tj^r 1 = 1 



Pm+2 m+2 

E n«:F'..4:'S 

^+27^?•m+2 = l U=l I 













't 




'''' ,w 

H-1 ■ ■ ■ 


■ ■ ■ 'ri+l 





• (2) 

h2 -(2) 

*g2-l ■ ■ ■ V2+I 



• (2) 

• (2) ^r2 



•(m+2) 

^5m+2 •(m+2) 

9m+2-l ■ ■ ■ '■rTO+2+1 



• (m+2) 
• (m+2) ^rm+2 



XT(X.ii) ...X.m ■■■X.(2) ...X. 



(2) 



'92 + 1 



{ m +2) ...X.( m+2} 
9m+2 + l ''m+2- 



1 



-^.((l®r)(rr[(//)-+2])), 



m + 2 

where in the next-to-last equality we note that when the operator ^ is applied, the various 
"boxes" are cyclically rearranged producing m + 2 terms which are identical after permuting 
the indices in the sums over pi, . . . ,Pm+2 and qi, . . . , qm+2- n 

Lemma 3.5. Assume that f = Qig, g = g* e and assume that 1 (g) 1 e domdj for all 

j. Assume that \\^ f\\A(s>^A < 1 and let 



Qig) 



:i ® r + r ® 1) 5^ (^/),, - (log(l + ^/)),, 



Then 



! + <// 



Proof. We use Lemma 3.4 and compare the (converging) Taylor series expansions of both 
sides term by term. □ 

Lemma 3.6. Assume that r is the semicircle law. Let 

K{f) = -/\/f)-f. 

Assume that f — QJg, and that g ^ is cyclically symmetric. Then 



K{f) 



|(l®r + r®l) {Y.[/f])j 



= 9{{l®T + T®l){Tr[/9g\)-J^g}. 
Proof. When m = —1, the equality in Lemma 3.4 becomes 

^ {(1 ® r + r ® 1) (Tr [/ f])} = -/* {/ f) + /f#/* (/) , 
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where I G Mnxn{-^ ® denotes the identity matrix. Since r is the semicircle law, 

= X, and ^ f^X = ,yV f (here ,yV is applied entrywise to the vector /). 
On the other hand, since ^ reduces the degree of polynomials by 1, JV = ^^g — ^g = 
^jVg — f. Thus we deduce 

^{{l0T + T(^l){Tr[^^g])-jKg} = - f) + ^ ^g - g 

= K{f). 

Using the equality /#X = ^g^X = ,yV g^ we get the statement of the Lemma. □ 

Lemma 3.7. Let r he the semicircle law. Assume that f = ^g for some g = g* cyclically 
symmetric and that \\^ f\\A<»„A < 1- Let 

Q{g) = [(1 ® r + r ® l)Tr {{^ ^g) - (log(l + ^ ^g))}] . 

Then Equation (3.4) is equivalent to the equation 

^{{l®T + T®l){Tr[J^g])- yKg} = ^{W{X + ^g)) + ^Q{g) + /^gi^^g. 

Proof. By Lemma 3.6, the left-hand side is precisely K{f). By Lemma 3.3, we obtain 

K{f) = nW{X + /)) + (//#^* il+J]) ~ ^* (l^^T?) ^ 
which proves the statement. □ 

We now turn to the proof of existence oi & g that satisfies the equation above. We will 
use a fixed-point argument; thus we first give some estimates on differential operators that 
will be involved. 

3.5. Technical estimates on certain differential operators. For the remainder of the 
section, we will assume that r : ^ ^ C satisfies 

(3.6) \r{q)\<Ct''' 

for any monomial q € . We begin with a few technical estimates on certain differential 

operators. 

Lemma 3.8. Let gi, ■ ■ ■ , gm ^ M)- Set 

Qm{9l, . . . , 5m) = (1 ® T + T ® 1) I ^ [{^ ^9l) ' ' ' {/^9m)\ 

I i=l 

Assume that (3.6) holds and moreover that Cq/A < 1/2. Then 

m 

||g„(E5i, . . . , T.gm)\\A < 2{2A-^r U IMa. 



k=l 



In particular, Qm extends to a bounded multilinear operator on s^^-^^ with values in s^^-^^ . 
Proof. Note that by definition, assuming that gf is a monomial, we have 

<-ieg g 

g=AXjB BA=RXiQ 



FREE MONOTONE TRANSPORT 21 

Thus ii gi, . . . , are monomials, then 

lit S^* BkAk=RkXj Qk 

and so, using the fact that for a given there can be at most degg^ decompositions of g^ 
as AkXji^^^Bk, and the degree of determines such decompositions, we conclude that each 
such sum has at most Ylt Qt nonzero terms, so that 

11(1 ® T)[{/&T.g,) ■ ■ ■ {J&T.g^)],,aA 



9k=^k^j^. Bk 



-2+- •+deg i/77 



<; J^i-\ Hm . (jdeg gi-h- 

/ \ deggi-ii— 2H |-deggm-im-2 

_ ^deggi+-+deg5,n^-2m ^ ( ) 



m ^ m 



k=i k=i 

A similar estimate holds for (r (g) ^T^gi) ■ ■ ■ ^^gm)]jojm- 

Since Qm{^gi, ■ ■ ■ , ^gm) is multi-linear in gi, . . . , gm, we find that if gj = ^ A(j, g)g is the 
decomposition of gj in terms of monomials, then we get that 

Qm{^gi,---,^gm) = 5Z K^,Qi)---K^^Qm)Qm{^qi,---,^qm) 

and so 



||g^(E^i,...,E^j|U < 2(2^-2)- E ni^(^'^^-)llk.-IU 

gi,—,gm 3=1 

m 

m 

= 2(2^-2)- nii^.iu- 

This concludes the proof for arbitrary gi, . . . ,gm- □ 

Lemma 3.9. For g, f E s^o, set Qjn{^g) — Qm{^9, ■ ■ ■ , ^g)- Assume that (3.6) holds and 
moreover that Cq/A < 1/2. Then 

m—l 

WQmiT^g) - gm(s/)iu < 2{2A-'r E ii^ii^ii/iir'"i/ - ^lu 



fe=0 



InpaHicular, ||Q„(E^)|U < 2{2A-^)"'\\g\\^. 
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Proof. Let k = 2{2A-^)"'. Then 



\\Qm{m-Qm{^9)\\A 



m—l 



m— fc k+1 m—k—1 



k=0 



m—l 



fe=0 

m—l 



m— fc— 1 



< '^Eii^II^II/IIa"'"'II/-^IU- 

fe=0 

□ 

CoroUciry 3.10. Assume that (3.6) holds and moreover that Cq/A < 1/2. Then the maps 
9 '->■ Qmi^g) eoctend by continuity to the completion of with respect to \\ ■ \\a- 

Lemma 3.11. Assume that (3.6) holds and moreover that Cq/A < 1/2. Let g G -^d^^ be 
such that \\g\\A < 4"? 



(—lY 

QiT.g) = Y,y^Qra^2{^g). 



m>0 



m + 2 



Then this series converges in \\ ■ \\a- Moreover, in the sense of analytic functional calculus 
on Mnxn{W*{£/ O ^"P, T O r°P)), we have equality 



Qi^g) 



;i ® T + T ® 1) i ^ U^J^g),, - (log(l + ^^s^)), 



i=i 



Furthermore, the function Q satisfies the (local) Lipschitz condition on {g : \\g\\A < A'^/2} 



||g(E^)-g(E/)|u<||/-5lU 



A' \ (i-^) (i-^) 



and the bound 



\\Qim\\A< 



2||g|U 



1 - 



2MU 

A2 



Proof Let k = and A = By Lemma 3.8, we have that ||Qr«+2(^)|U < 2(A//t)('"+2) 

. Thus, as we assumed A < k, the series defining Q(llg) converges. To see the claimed 
equahty with the expression given by functional calculus, we only need to note that log(l + 
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x) — — J2m>i ■ Finally, since m + 2 > 2 in our series, we obtain that 

||g(s^)-g(E/)iu < 5^^-||g^+2(s^)-gm+2(s/)|u 



m+l 



m>0 fe=0 

< 11/ - 9U2A-') ( E E(2^"')'ll/ll^(2^"')'ll^ll^ - ' 

\l>0 k>0 / 



where we have written m = I + k — 1 which is non-negative precisely when I and k are not 
both zero. Thus, if ll/IU < and II^IU < we deduce that 

WQim - QmU < 11/ - ^lU^ I 7 — ^^^^7 — - 1 I • 



Setting / = gives us the estimate 



lig(Ej)iu < ?!lf|^ (^^-^ - 1 1 , 

as claimed. □ 

CoroUairy 3.12. Let g e assume that \\g\\A < and let B > A + \\g\\A- Let 

W e =2/^^). Assume that (3.6) holds and moreover that Cq/A < 1/2. Let 

Fig) = -W{X)+ (^{l^T + T^l){Tr{^S'Eg)) 

- i^W{X + 9T.g) - W{X) + Q{Y.g) + ^^E^#^E^| 

= -WiX + ^Hg) - ^^Ey#^Ey + (l(8)T + T® l)rrlog(l + ^^E^). 

(Tiere u;e abbreviate Vr(Xi, . . . 6?/ etc.j 

r/ien F[g) is a well-defined function from =2^^"^^ to Moreover g i->- F(g) is locally 

Lipschitz on {g : \\g\\A < A^/2}: 

iiFfe) - F(/)iu < 11/ - ,iu{:^ ( + 1 



A2 



.) (l - 3^) 



11 

+ \\JW\U^B + -{\\f\\A+\\9\\A) 



and bounded: 



\\n9)\\A < \\g\\A\ ^ (^YT^ + ^ 1 + W/mU.B + ?||5lU ^ + WWa 
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In particular, if 
(3.7) 

then F takes the hall {g 
X < I on that ball. 



A> A, <p<l 



\\/W\\ 



A <^ -P- 

A ^ 12n 



{A+p)®-„{A+p) 



< 



qWa < £'} itself and is uniformly contractive with constant 



Proof. Note that ||^S5f||^ < \\g\\A and so ||X + ^S^fH < B. The Lipschitz property then 
follows from Lemmas 3.11 and 3.9, as well as the estimates: 

\\W{X + ^T.g)-W{X + ^^f)\\A < \\^W\\B(,^B\m9-^m\A 

< \\^W\UM\g-f\\A 



and 



< -\mg4fmg - ^Ef)\\A 



--||^E/#(^E^-^E/)|U 



n , 



< 



n 

^\\9 



fhmiA 



The estimate on ||F((yf)||^ follows from the identity -F(O) = —W and the Lipschitz property. 
Assuming that (3.7) holds, so that A > 4 and ||gf||yi<^<p<l, we have that 




1 - 



yl2 



1 - 



A2 



+ 1 < 



64 \ 1 
49 + ^ <2 



and so the Lipschitz constant of F is bounded by ^+n^ + \\^W\\(^A+p)(s>^{A+p) < 
Also, 



1 + 1. 



< 



P 1 



4n 8 V 7 



1 1 



\W\\a< 



P_ 
An 



P_ 
An 



so that the image of the ball {g : \\g\\A < £-} is contained in that ball. 



□ 



3.6. Existence of g. We will now consider the case that r is the semicircle law. Then (3.6) 
holds with Co = 2. We remind the reader that H, E, are defined in section 2.2 whereas is 
defined in Definition 2.2. 



(A) 



IS 



Proposition 3.13. Let r be the semicircle law. Assume that for some A, W e < 
cyclically symmetric and that conditions (3.7) are satisfied, i.e., 

A>A, 0< p<l 

\mA<^ 

\\^W\\^A+p)^^{A+p) < I 

Then there exists g and g = T,g (so that g = ,yVg ) with the following properties: 

(i) Both g and g belong to the completion of with respect to the norm \\ ■ \\a; 

(ii) g satisfies the equation g = ,S^Y{F{g) 

(Hi) g and g depend analytically on W , in the following sense: if the maps (3 h- )■ Wp are 
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W{X + 9g) - l^g#^g + (1 r + r ® l)rrlog(l + 



analytic, then also the maps (3 i->- g{f3) and (3 ^ g are analytic, and g ^ if \\W\\a — >■ 
(iv) g satisfies the equation 

(3.8) yU^g = -W{X) + yu(^{l ^ T + T ^ l){Tr{^ ^g))- 

i^W{X + %) - WiX) + Q{g) + \^g4^g 

or, equivalently, 

(3.9) - yU^g + ®T + T® l){Tr{J^g)) = yiiW{X) 

+ |l^(X + %) - W{X) + Q{g) + ^%#% 

Proof. Observe that Equation (3.8) is equivalent to 

yu^g = yUF{^g) 

with F the function defined in Corollary 3.12, where we used that W is cychcally symmetric 
and in ^q^^ Existence of g follows, essentially, from the Imphcit function Theorem 2.3. We 
prefer to repeat its proof here for completeness. We set = W{Xi, . . . , X„), and for each 
A; > 0, 

gk = yUF{gk^,), 

Since is a linear contraction, the last part of Corollary 3.12 implies that under our 
hypothesis, ^11 o F is uniformly contractive with constant 7/8 on the ball B — {g : \\g\\A < 
^}, and takes this ball to itself. It follows that gk ^ B for all k, and moreover 

7 

\\gk - gk-i\\A = \\yn.F{gk-i) - yn.F{gk-2)\\ < -\\gk-2 - gksWA, 

so that gkS converge in || • \\a to a fixed point g. Since g e 'S^o^^ , so does g — "Eg. This 
proves (i) and (ii). 

Since (assuming that /3 i— )> IV^ is analytic) each iterate c/k is clearly analytic in /3 and the 

convergence of the Banach-space valued function gk{/3) G ^q^^ is uniform on any compact 
disk inside \f3\ < (3o, it follows from the Cauchy integral formula that the limit is also an 
analytic function. It follows from the norm estimates in Corollary 3.12 (taking p — )■ 0) that 
the bound on ||^fe|U iterative procedure goes to zero as — >■ 0, so that \\g\\A — >■ 

as \\W\\a 0. 

Part (iv) follows from the definition of F and S. □ 

We leave the following to the reader (note that g is clearly cyclically symmetric, and it is 
not hard to see that this implies that also g is cyclically symmetric): 

Proposition 3.14. Assume that W = W* and that W satisfies the hypothesis of Proposi- 
tion 3.13. Let g be the solution to (3.9) constructed in Proposition 3.13. Then g belongs to 
the closure of cyclically symmetric polynomials and also satisfies g = g*. 
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3.7. The map / = Sig satisfies Equation (3.4). Prom now on, we will assume that r is 
the semicircle law. 

Theorem 3.15. Let A> A' > A. Then there exists a constant C{A^A') > depending only 
on A, A' and n, so that whenever W = W* e s^^^^^^ satisfies \\W\\a+i < C{A,A'), there 
exists an f E {s^^"^ which satisfies (3.4)- In addition, f = 2ig for some g G s^'^'^ \ The 
solution f — fw satisfies \\fw\\A' — >■ as — >■ 0. Moreover, ifWp is a family which 
is analytic in (3 then also the solutions fwp are analytic in (3. 

Proof. We first note that, by the Cauchy integral formula, if we choose Ai so that A < Ai < A, 
then 

M n(A + 1) „ „ 

Similarly, we see that 

(3.10) IU/IU,«.A, < ^^^^^^II^IU 

and also 

\\^w\\^A,+mAM+i) < ^^^ll^ll^+i- 

Thus we can choose C{A,A') > so that if < C{A,A'), then the hypothesis of 

Proposition 3.13 is satisfied, and there exists some g which satisfies 

- yU^g + yU{l ® r + r ® l)(Tr(^%)) 

= yUW{X) + yu !^W{X + %) - W{X) + Q{g) + ^%#%| . 

Purthermore, we may assume (by choosing a perhaps smaller C{A, A')) that / = ^g satisfies 
ll<?^/IU'®^A' < 1- Applying ^ to both sides and noting that for any h, ^yUh — Sih, gives 
us 

9{{l®T + T®l){Tr[J9g\)-^g} = ^{W{X + f)) + 9Q{g) + J f4f, 

which, according to Lemma 3.7 is equivalent to (3.4). The hypothesis of Lemma 3.7 is 
satisfied since the map T : ^ i— )■ ^ fj^i, is strictly contractive on {^s^^^ "^')^ because of the 
bound on ^ f (and so the map id+T has no kernel). □ 

Since IV = X ^ fw — )■ X is bounded by A as — in it follows that for 

sufficiently small || the law of Y is the unique solution to the Schwinger-Dyson equation 
with self-adjoint potential V = ^^Xj + W, and thus the law of Y is exactly Ty (by 
[GMS06, Theorem 2.1] or alternatively since, for sufficiently small, V is then (c, M) 

convex and the Schwinger-Dyson equation has a unique solution [GS09]). In particular, 
W^*(Yi, . . . , Yn) = W*{rv)- We thus obtain a free analog of Brenier's theorem: 

Theorem 3.16. Let A > A' > 4, and let Xi, . . . , X^ G (M.r) be semicircular variables. 
Then there exists a universal constant C = C{n, A, A') > so that whenever W € s^^"^^^^ 
satisfies < C, there is an element G E s^^^ ^ so that 

(Fi,...,r„) = (^iG',...,^„G')e^/(^') 

has law Ty,V = \Y, X] -f- W. 

Moreover, the Hessian ^ Q)g is a strictly positive element of Mnxn{M^M°P). 
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In particular, there are trace-preserving injections C*{tv) C C*g^(F„) and W*{tv) C 

If the map /3 ^ Wp is analytic, then Yi,...,y„ are also analytic in (3. Furthermore, 
\\Yj — XjWa' vanishes as \\W\\a+i goes to zero. 

Proof. We use Theorem 3.15 and set G = | ^ X| + (7 and the discussion before the statement 
of the present theorem. The only thing left to prove is that the Hessian of G is strictly 
positive. But the Hessian of G is given by 

J^G = 1 + J^g. 

In the proof of Theorem 3.15, we have chosen C in such a way that S!g^A'®.^A' < li 
which means that ^ Q)G is strictly positive. □ 

4. Some applications. 

4.1. Analyticity of the solution to the Schwinger-Dyson equation. Under the hy- 
pothesis of Theorem 3.16, we deduce that if /3 — ?> is an analytic family of potentials and 

= I ^ X| + Wfi^ then there exists an analytic family of elements Yg whose law rg satisfies 
the Schwinger-Dyson equation for Vg . We deduce the following corollary, which was already 
proved in [GS09], by a different method: 

CoroUciry 4.1. Let P ^ he a fixed polynomial. Then Tp{P^ is analytic in fi in a neigh- 
borhood of the origin. 

4.2. Isomorphism results. Applying the implicit function theorem, we can improve The- 
orem 3.16 at the expense of possibly choosing a smaller bound on W . Indeed, we let 
(Yi, . . . , Yn) = f = ^g be as in Theorem 3.16, we can always assume that y,(0, . . . , 0) = 
by replacing each Yj with Yj — Yj{0, . . . ,0). Thus we may assume that /(O) = 0. 

Theorem 4.2. Let A > A' > 4, and let Xi, . . . , Xn G (M, r) be semicircular variables. 
Then there exists a universal constant C — C{n, A, A') > so that whenever W e ^/^'^'^^^ 
satisfies \\W\\a+i < C, there is an element G e j^/^^'^ so that: 

(1) // we set Yj = ^jG, then Yi, . . . , y„ G has law Ty, with V^lJ2Xj+ W; 

(2) Xj = Hj{Yi,..., Yn) for some H G 

(3) the Hessian ^ is a strictly positive element of Mnxn{M®M°^). 
In particular, there are trace-preserving isomorphisms 

C*(ry) = C*(Xi,...,X„), W*{Ty)=L{¥n). 

4.2.1. Isomorphism class of algebras with analytic conjugate variables that are close to the 
generators. 

Corollary 4.3. Let Xi, . . . ,X„ G (M, r) be generators of a von Neumann algebra M and 
assume that the Fisher information . . . is finite. Assume moreover that ^* — 

Cj = d*{Xj) belongs to for some A > A. 

Then there exists a universal constant C = C{n,A) so that if — XjWa+i < C, then 
M = L(Wn) and also C*{Xi, . . . , Xn) — C*{Si, . . . , Sn) where Si, . . . ,Sn are free semicircular 
variables. 

Proof Letting V = (Ei=i -^iO + ' one gets that Q = (cf. [VoiOO]). Thus 

we can write V = ^Y^Xj +W; then for small enough \\Q — Xj\\A+i, W = W* satisfies the 
hypothesis of Theorem 4.2. □ 
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4.2.2. Voiculescu's conjecture with polynomial potentials. The following corollary is a partial 
answer to a conjecture of Voiculescu [Voi06, p. 240] (the full statement of Voiculescu's 
conjecture is that the isomorphism should hold for arbitrary polynomials W and arbitrary 
/3; however, even uniqueness of a free Gibbs state is unknown in that generality!) 

Corollary 4.4. Let A > A and let W e ^^"^^ be a self-adjoint power series. Let Vp = 
J2 + Oj'f^d let Tff be a trace satisfying the Schwinger- Dyson equation with potential 
Vp. Then there exists a > Q so that W*{Tp) = Vr*(ro) = L(F„) and C*{Tp) ^ C*(to) = 
C*(S'i, . . . , Sn) for all —f^o < (3 < (here Si, . . . ,Sn is a free semicircular system). 

4.2.3. q-deformed free group factors. In [BS91], Bozejko and Speicher have introduced a 
family of von Neumann algebras ry(]R"), which are "g-dcformations" of free group factors. 
For q = 0, rq(IR.") = L(F„), but the question of whether an isomorphism like this holds 
for g 7^ remained open. Despite being much-studied ([BS91, BKS97, Nou04, SniOl, Sni04, 
Shl04, Ric05, Shl09, KNll, DablO, Avsll] is an incomplete hst of results about these factors), 
the question of the isomorphism class of these algebras for all values of q remains elusive. 
Nonetheless, we are able to settle it for small q (depending on n < oo): 

Corollary 4.5. Let TglMT') be the von Neumann algebra generated by a q- semicircular system 
S^i \ . . . , Sn \ n G {2, 3, . . . }. Then there exists a < go = Qoin) depending on n, so that 
r,(M") = ro(M") = L{¥n) and C*{S[''\ si'^) = C*{S[°\ ) for all \q\ < go- 

Proof. This follows from the fact [DablO, Theorem 34] that the conjugate variables = 
q = to {S{'^\ SL^'^) exist and that for some A > 5, - XjWa ^ as g ^ 0. □ 

4.2.4. Free entropy. Note that under the assumptions of Theorem 4.2, we end up expressing 
Xj's as a convergent power series in n semicircular variables Si, . . . , Sn- By the change of 
variables formula of Voiculescu [Voi94], we get: 

Corollary 4.6. Let V = \ ^X'j + W , and assume the hypothesis of Theorem 4-2. Let Ty 
be the free Gibbs state with potential V. Then 

X{rv) = X{r) + r ® T[Tr(log(^^G))] 

where r is the semicircle law and G is as in Theorem 4-2- 

4.3. Monotone transport for random matrix models. Let V = \ ^X| + W^ and Ty be 

the free Gibbs state with potential V . Assume that W satisfies the hypothesis of Theorem 
4.2 for some A> A' > A and let F = (Fi, . . . , F„) = &G be the map in (^(^'))" constructed 
in that theorem. We let 

l^(Xi,...,X„) = i^X| + W^(Xi,--- ,X„). 

on {maxj < A} and V = +oo otherwise. V is strictly convex for small enough 

Let be the measure on (M^^jy)" given by 

d/^if ^ = ^ eM-NTr{V{Ai, Ar,)))dAi ...dA^ 

where we set dAi . . . dA^ = Lebesgue measure. Consider also the Gaussian measure 

V^) = ^exp 



-NTr 



dAi . . . dAr 
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We denote by x : R — )■ [—A', A'] a smooth cutoff function so that x(^) = for a; G [—4, 4], 
|x(a;)| < A' and x has a uniformly bounded derivative. On the space of n-tuples of N x N 
self-adjoint matrices (Af^Txiv)" with norm no more than 4, consider the functional calculus 
map 

■.(A,,...,Ar.)^ (FjixiA,), x{An)))]=,. 
Finally, let F^^^ — V^'^^^ be the unique monotone transport map on (Mf^^j^)'^ so that 

Theorem 4.7. Let A > A' > 4 and assume that \\W\\a+i < c with c small enough, at least 
smaller than the constant of Theorem 4-2- 
(i) With the above notation, 



hm [ -Tr ( [f(^) - F^^A d^i^''^ = 0. 



In other words, as N ^ oo, the "entrywise" monotone transport F^^^ is well- approximated 

as N ^ oo by the "matricial" functional calculus map F^^^ . 

(a) Assume in addition that W is {c, A) -convex for some c > 0. Then ^ F < 1. Further- 
more, j^Tr ^ F(^) — F(^) ^ vanishes almost surely as N goes to infinity. 

Proof. Let us put, for a probability measure v on with density p(a;i, . . . , XM)dxi . . . dxM, 
and U : ^RU{+oo} a convex function, 

H{u) = j p{xi, xm) \ogp{xi, XM)dxi ■ ■ ■ dxM 

Hu{iy) = J p{xi,...,XM)'fogp{xi,...,XM)dxi...dxM + J Udu. 
If h is an invertible transformation with positive-definite Jacobian Jac h, then 

= -ff(z^) — y logdet Jac/i dv 

= H{iy) - J Trlog Jac/i du. 

Both entropies Hu[v) and H^v) are convex functions of the density of v. Moreover, Hu{h^i>) 
and H{h*v) are convex vahiiU is convex. 

It is well-known that Hu{y) is minimized precisely by the Gibbs measure with potential 

U: 

Huiv) ^ vcd Huiv') <^=^ dv{xi,. . . ,Xm) ^ -;^Qy^Y>{-U{xi,. . . ,XM))dxi. . .dxM- 
Similarly, for V € let x(r) be Voiculescu's microstates free entropy, and put Xvl^) = 

x{j)-T{y). 

Then using [Voi02b, GS09, Gui09] one has 

Xvij) — supxv('7'') <^=> T is the free Gibbs law with potential V. 

t' 

Moreover, if we set V^^\Ai, . . . , A„) = NTr{V{Ai, . . . , Ajv)), then following [GS09, proof 
of Theorem 5.1] or [GMS06, Section 3.3] we see that /x[^\;^Tr(P)) converges towards tv{P) 
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for all polynomials P, in particular the limit does not depend on the cutoff provided ||iy 
is small enough. In fact, has a strictly log-concave density, Brascamp-Lieb inequal- 

ities allow one to show that the matrices under //^■* are bounded by 4 with overwhelming 
probability. Moreover 

1 77/ 

Xv{rv) = limJ—\ogZ^ + -\ogN) 



^im(|logiV-i,i/^(.)(;.r))) 



Furthermore, we claim that, if jly — Fi^^//^^\we have 

1 

To see this, we first note that 



lim -^Hy(N){fl'^^] 



-Xv{rv). 



Moreover, 



1 



JV) 



-^Tr log JacF(^) 



On the other hand (see e.g. [GuiOQ]), because of concentration phenomena, as //^"^^ has a 
strictly log-concave density, if we take r to be the semicircle law, then 



— rTrlogJacF^^) 



^ T®T(Tr log ^F). 



Thus 



Hijj'y ^) ^ -x(t) - r ^ T{Tr log ^ F) = -x{tv) 



so that using Corollary 4.6, 



ji^^^yw(/iif^) = -Xv{rv) 



as claimed. By convexity of entropy as a function of the transport maps, for e e [0, 1], 



< (1 - 5)i/y,.v)(Ff V^"^^) + £i?yw(Fi^V(^)) - Hyi^^ifl 

= e{Hy(N) (/iSf )) - Hy(N) (/i^f ^ ) ) . 

On the other hand, let A^^) = F(^) - F(^) and set 



1/ ; 



Since i/v-c^) ( [(1 - ^)F^''^ + ^F^""^] P^^""^) 

has an absolute minimum at £ = 0, we deduce 

that Dtv > 0. 

Furthermore, assuming that is small enough so that on the set {maxi<j<„ II Alloc ^ 

^4} the Hessian Hess(V*^^)) is bounded from below by cN times the identity operator for some 
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c > 0, and noting that by definition the image of F^^^ is bounded by A whereas the image 
of F^^^ composed with x is also bounded by A, 

Hy,., ( [(1 - £)F W + £F ^ - sD^ - Hy,., (/.if)) 



(4.1) 



= E (jv) 



E„(jv) 



tl Tr 



(jac(iFW + (l-t)FW) 



Jac(AW) 



+rr(Hess(yW(FW +iAW))(A^)2) Li 



> 



^c7VE^(.)(Tr((AW)^). 



We thus get (recahing that Djq > 0) 

< £(//vw(/iif))-//^w(//if))). 



Since 



HyiN){fi\^^) - Hy(N){lj!'y^) 



we can now choose e — sJoln — > to conclude that 



^ 



(4.2) 



N 



< a 



1/2 
N 



which completes the proof of the first point. 

Let us now assume that W is (c, 74)-convex for some c > 0. Let l^(^) — NTr{W{Ai, . . ., A„))| 
on matrices satisfying maxj IjAjHoo < A' and infinite otherwise. Then W^^^ is convex. It 
then follows from Caffarelli's results [CafOO, Vil03] that the optimal transport map F^^'' 
taking the Gaussian measure n^^^ to the measure with density Z]^^ exp{—^N J2Tr{A'j) + 
W'^^^Ai, . . . , ^4^)) has Jacobian uniformly bounded by 1. 

The map (^i, . . . , An) — > j^Tr{{A'^^\Ai, . . . , An)"^) can be viewed as the composition of 
the map R:{Ai,..., aJ) AW(Ai, . . . , A„) and the map Q : {Ai, . . . , A„) ^ ^ ^ Tr{A]). 
The Jacobian of R is given by - Hess^(^), where ^(^^ is such that = V^^^). 



If maXi 



l^jlloo < A, 



this is bounded (as an operator on Hilbert spaces 



NxN 



— ^ M" 



endowed Hilbert Schmidt norms \\A\\ = ^jTr{A*Aj)) because is a power series and 
because of Caffarelli's bound. Hence the map R is Lipschitz with a uniform Lipschitz constant 
on the set where maxj ||Aj||oo < A. The map Q is Lipschitz with Lipschitz constant of the 
form C/VN (see [Gui09, Lemma 6.2]). 

Therefore, by concentration inequalities and (4.2), we deduce 



lim lTr((AW)') = // 



(JV)_ 



a.s. 



We next come back to (4.1) and observe that t Jac(-F*-^-') is bounded above uniformly (as 
a small smooth perturbation of the identity) as well as (1 — t) Jac(F^^'') which is bounded 
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above by one by [CafOO], so that we get for £ < 1/2 



Taking once again e — s/on — > 0, we obtain that 



1^' 



^Tr I (jac(F(^))-Jac(F(^)))' 



^ 0. 



We now apply [CafOO] to conclude that 
(4.3) JacF(^) < 1. 

Let Yn — Jac(F(-^)), = Jac(F(^)) be random variables taking values in the space 
M„xn(End(M^"^jv)) — MnxniMNxN<^MNxN) cudowcd with the normalized trace -^Tr^Tr. 
Both Yjsf and Y/v are bounded in operator norm, and consequently define elements Y and Y in 
the ultraproduct von Neumann algebra Mnxni^NxN ® Mnxn)- Since ||liv — ^Ar||2 0, 
Y = Y. 

Once again, because of concentration. 



i^Tr{(Jac(F(^)))^} 



^ T (g) TTr{{^F)P) 



so that the spectrum of ^ F is the same as that of Y (and so the same as that of Y). 

By (4.3), the spectrum of Y^ lies in the interval [0,1]. But this imphes that also the 
spectrum of the limiting operator Y is contained in the same set. Thus ^ F <\. □ 

5. Open questions. 

We list some open questions that are raised by our results. 

(1) In the classical case, Brenier's theorem [Bre91, Vil03] asserts much more than the 
statement of our main theorem: the classical analog of the map 3)g gives optimal 
transport from r to Ty for quadratic Wasserstein distance. It would be nice to under- 
stand if the same holds true in the non- commutative case (see [BVOl] for the exten- 
sion of the notion of the Wasserstein distance to non- commutative random variables). 
Note that the map we construct is optimal in the single- variable case n = 1. 

(2) Brenier gave a heuristic derivation of his theorem through a very general "polar factor- 
ization" theorem. Does a theorem like that hold in the non-commutative case? There 
is an infinitesimal analog of his decomposition (related to the classical Helmholtz de- 
composition of vector fields) which has been extensvively studied by Voiculcscu in 
[Voi02a]. Can non-commutative monotone transport be also obtained in the same 
way? 

(3) Does the positivity condition ^ 9G e {F > : F € M„x„(r(8)T°P)} on the "Hessian" 
of G translate into any kind of convexity properties of Gl 

(4) What happens to our map F = 2iG in the case that V is not strictly convex? 
It can be seen that in the absence of the bounds on W the isomorphism (or even 
the embedding C*{tv) C C*(5'i, . . . , Sn)) fails to exist at least on C*-level. This is 
due to the fact that for certain non-convex polynomials, solutions to the Schwinger- 
Dyson equations may lead to C*-algebras with non-trivial projections. But the free 
semicircular system generates a projectionless C*-algebra. Thus failure of convexity 
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oiV~ ^^Xj +W must be 'Visible" as a defect of regularity of the transport map 
F. 

(5) Is the monotone transport map unique? More precisely, let Xi, . . . , X„ be a semicir- 
cular family, and assume that Y — (1^, . . . , y„) and Y' — (Y/, . . . , Y^) both belong to 
the L^-closure of {^g{Xi, . . . , Xn) : > 0}. If the law of Y is the same as the 

law of r' , is y = Y'l 

Wc can prove a uniqueness statement for our monotone transport if we assume more on the 
transport map and the "target" n-tuple Yi, . . . , F^: 

Theorem 5.1. Let Xi, . . . , X„ be a semicircular family, and assume that Y — (Yi, . . . , Y^) 
and Y' = (Y-[, . . . , Y^[) both belong to the L'^-closure of {^g : ^g > 0}. Assume moreover 
that Y and Y' are both invertible non- commutative power series in Xi, . . . , Xn- IfY and Y' 
have as their laws the same free Gibbs law Ty with (c, M) -convex potential V , then Y = Y' . 

Proof. The assumption that Y, Y' are invertible non-commutative power series allow us to 
apply Voiculescu's change of variable formula for free entropy [Voi94]: 

X{Y) = x{X)+T®T{Tr{\og/Y)) 
X{Y') = x{X)+T®T{Tr{\og/Y')). 

Since Y and Y' have the same law, they have the same free entropy (which is finite [GS09]). 
Thus we conclude that 

T ® r(Tr(log JY)) - t{V{Y)) = t® r(Tr(log /Y')) - riViY')). 

Let 

i:{Y)^x{Y)-x{X)-T{V{Y)). 

Then is maximal iff Y has the law ry (see e.g. [Voi02b], we sketch the argument for 
completeness: if one replaces F by F -|- eP{Y) for some polynomials Pi, . . . , P„, then 

V'(F + eP{Y)) = z^(y) + e {r ® T{Tr{/YP)) - t{9V{Y)P)} + O(e'), 

and so any maximizer to ^ satisfies the Schwinger-Dyson equation and (because V is c, M- 
convex), has the same law as F.) 
It follows that 

r ® r(Tr(log /Y)) - t{V{Y)) = x(F) - r(V(F)) = maxx(^) - r(V(Z)), 
so a fortiori 



T (g) T(rr(log ^F)) - t{V{Y)) = sup r (8) T(rr(log jfY)) - t{V{Y)). 

Ye{Sig:J'&g>0} 

But since r (8) r(Tr(log ^Y)) — r(y(Y)) is strictly convex in F for ^Y in the positive 
cone of Mnxn{M®M"P) (here M = W*lXi, . . . ,X„)), it follows that there is at most one 
n-tuple F in the closure of F e {^g : ^ ^g > 0} which gives this maximal value. Thus 
F = F'. □ 
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